<html>
<head><title>����</title></head>
<body>
<center>
    <table width=780 border=0 cellpadding=3 cellspacing=0>
        <tr>
            <td style='line-height:150%;font-size:16px;text-indent: 20'>
                <p><a href="readme.htm">����</a></p>

                <p align=center><b>�ھ��� ��ʳԭ��</b></p>

                <p><b>һ����Ӱ�ᡢ��Ӱ��α��Ӱ����Ӱ����ʳ����</b></p>

                <p>
                    ������ס��̫���⣬ʹ�����ڵ����ǿ�����̫����ֻ����̫����һ���֣����γ�����ʳ��������ʱ�̣�������Ӱ�����ڴ��ڵģ��������������Ӱ�����򣬲��ܿ�����ʳ����</p>
                <center><img src="img/a012.gif"></center>
                <p>
                    ���Ƕ�����һ������ʱ���ܷ�����ʳ��������Ӧ�ڵ����̫��֮��ſ��ܷ�����ʳ�������º�˷ʱ�̣���˿̸���ſ��ܷ�����ʳ��������º�˷ʱ������λ��̫���ߡ���̫���͡���Ӱ�����ڵ�����Ϸ�������֮�ϣ����·����ϼ�֮�£��������Ͽ�������ʳ��</p>

                <p>
                    ����������һ������ʱ���ܷ���ȫʳ�������ĳ���뾶Ϊ1����ô������Ӱ�ӳ���ԼΪ58.57�������������Χ����ңԶ��̫����ת����ת�Ĺ���Ǹ���Բ��������������ʱ��̫����ˣ���ʱ��ԶһЩ������Ӱ�ӳ���ԼΪ58.57��1.67%�����ֵ����¾���60.34�ܽӽ������Ƶ��˶�����Ǹ�ƫ���ʽϴ����Բ�����¾�60.34��5.5%����˵��¾������С��Ӱ�ӳ��ȣ��γ���ʳ��</p>

                <p><b>������ʳ��������</b></p>


                <p><b>��ռ�ֱ�������Ľ���</b></p>

                <p>
                    �������״�ɽ���Ϊ��ת������Ӱ�����Ӱ׶��ĸ�߿���ʹ��һ��ֱ����������ĳʱ����Ӱ�������Ľ�����Ǹ�ʱ����ʳ�����ĵ㡣Ӱ׶�ĸ���ĸ�������Ľ��㹹��Ӱ���ڵ����ϵ������ߡ�</p>
                <center><img src="img/a013.gif"></center>
                <p>����ֱ�ߵĲ�������Ƚ�����</p>

                <center><img src="img/a009.gif"></center>
                <p>����ʾ��������ֱ�������Ľ��㣬ֱ�߾���(x1,y1,z1)��(x2,y2,z2)���㣬eΪ���ྶ�ȣ�rΪ��뾶��<br>
                    function lineEll(x1,y1,z1, x2,y2,z2, e,r){ //��ռ��������������Ľ���(�����x1�Ľ���)<br>
                    &nbsp; var dx=x2-x1, dy=y2-y1, dz=z2-z1, e2=e*e, A,B,C,D,R,t, p=new Object();<br>
                    &nbsp; A = dx*dx + dy*dy + dz*dz/e2;<br>
                    &nbsp; B = x1*dx + y1*dy + z1*dz/e2;<br>
                    &nbsp; C = x1*x1 + y1*y1 + z1*z1/e2 - r*r;<br>
                    &nbsp; p.D = B*B-A*C; if(p.D&lt;0) return p; //�б�ʽС��0�޽�<br>
                    &nbsp; D = sqrt(p.D); if(B&lt;0) D = -D; //ֻ�󿿽�x1�Ľ���<br>
                    &nbsp; t = (-B+D)/A;<br>
                    &nbsp; p.x=x1+dx*t, p.y=y1+dy*t, p.z=z1+dz*t;<br>
                    &nbsp; R = sqrt(dx*dx + dy*dy + dz*dz);<br>
                    &nbsp; p.R1 = R*abs(t), p.R2 = R*abs(t-1); //��λ��ǧ��<br>
                    &nbsp; return p;<br>}
                </p>

                <p>�������(x,y,z)����󣬽���ľ���<font color=red>�� = atan2(y,x)</font>�����?��<font color=red> L = �� -
                    g</font>������γ�Ȧ�'Ϊ<font color=red>tan��' = z / sqrt(x*x+y*y)</font>������γ�Ȧ�'�����γ�ȦյĹ�ϵΪ��<font
                        color=red>tan��' = (b/a)<sup>2</sup>tan��</font>��ʽ��b/a�Ǽ��뾶�ȳ���뾶����Ȼ����γ��С�ڵ���γ�ȣ����ϱ����ͳ���϶�����ȡ�
                </p>

                <p>��󣬵��?��Ӧ������-180�㵽+180��,����γ��Ӧ������-90�㵽90�㡣���������У����?������������Ϊ����ô��L = g - ��</p>

                <p><b>�ġ���ʳ�������õ����ϵ</b></p>

                <p>��û��Ӳ�Թ涨Ҫ���������͵����ϵ��һ���˵����ȷ����ʳ�����У������õ�ֱ�����ϵ��</p>

                <p><b>��һ�����任</b><br>

                    ������תֱ�����llr2xyz()��<br>
                    function llr2xyz(JW){<br>
                    &nbsp; var r=new Array(),J=JW[0],W=JW[1],R=JW[2];<br>
                    &nbsp; r[0]=R*cos(W)*cos(J);<br>
                    &nbsp; r[1]=R*cos(W)*sin(J);<br>
                    &nbsp; r[2]=R*sin(W);<br>
                    &nbsp; return r;<br>
                    }<br>

                    ��ֱ��ת�������xyz2llr()<br>
                    function xyz2llr(xyz){<br>
                    &nbsp; var r=new Array(), x=xyz[0],y=xyz[1],z=xyz[2];<br>
                    &nbsp; r[2] = sqrt(x*x+y*y+z*z);<br>
                    &nbsp; r[1] = asin(z/r[2]);<br>
                    &nbsp; r[0] = rad2mrad( atan2(y,x) );<br>
                    &nbsp; return r;<br>
                    }<br>

                    �����������ת����llrConv()��<br>
                    function llrConv(A,E){<br>
                    &nbsp; var B=new Array(),J=A[0],W=A[1];<br>
                    &nbsp; B[0]=rad2mrad( atan2(sin(J)*cos(E) - tan(W)*sin(E), cos(J)) );<br>
                    &nbsp; B[1]=asin ( cos(E)*sin(W) + sin(E)*cos(W)*sin(J) );<br>
                    &nbsp; B[2]=A[2];<br>
                    &nbsp; return B;<br>
                    }<br>

                    ��˵����ת�������е���rad2mrad()����������ǽ��Ƕ�ת����0��2PI<br>
                    ���÷���B = llrConv(A,E)<br>
                    �����ܣ���A�����ת��B����꣬EΪA����B���ļнǡ�<br>
                    ��A��B�������ΪA��B�Ľ��㣬���������ͽ�����֮�֣���������û���ѡ����Ӧע�⣺<font color=red>��A�����ؾ��Ȳ��������ߣ����������ʱ�����γ�����ӣ����ʾ�����Ϊ��㣬��ʱEȡ��֮Ϊ����Eȡ����</font><br>
                    ����1��AΪ�Ƶ���꣬ת����BΪ�����꣬���������Ϊ����㣬��ô��A��(���Ƶ�)��������ʱ����γ�����ӵģ�����Eȡ��E�ǻƳཻ�ǡ�<br>
                    ����2��AΪ�����꣬ת����BΪ�Ƶ���꣬���������Ϊ����㣬��ô��A��(�����)��������ʱ����γ�����ӵģ�����Eȡ����E�ǻƳཻ�ǡ�
                </p>

                <p><b>��������������ϵ����任</b></p>

                <p><b>��������ϵ�Ķ��壺</b></p>

                <p>
                    z�ᾭ�������ƽ������Ӱ�ᣬ��̫���Ǹ�����Ϊ��ֱ����Ӱ���Ҿ�����ĵ�ƽ��Ϊ���档����������Ľ���Ϊx�ᣬxyz���ṹ���������ϵ�����涨x������ĳྭ��z��ྭ��90�ȶ���-90�ȣ�����y�������ڱ�����</p>
                <center><img src="img/a003.gif"></center>
                <p>��z����ĳһ�����B������дΪֱ�������ʽ����B[0,1,2]��jdʱ��̫��������ΪS[0,1,2]������������ΪM[0,1,2]�����ڹ�ϵʽB
                    = S -
                    M��Ӧע�⣬�����ϵʽ��һ��ʸ�������ϵ����S��MΪ�������ʱ�����Ƚ�����תΪֱ����꣬Ȼ�������㡣���õ���B������תΪ������ꡣ��ˣ����������ϵ�У�z�������ĳ������ʾ��<font
                            color=red>B = xyz2llr( llr2xyz(S) - llr2xyz(M) )</font>��ʵ���ϣ�B�ķ���ʮ�ֽӽ�̫��������������ʳ�����ڼ䣬B��S�ĽǶȲ��10���롣
                </p>

                <p><b>��������任��Ҫ�أ�</b></p>

                <p>���������������Ϳ����ڱ������ꡢ�����ꡢ�������֮��������任��<br>1�����ཻ�Ǧ�<br>2�����ྭ(A��ĳྭ,�����ྭ)��<br>3���������κ���ʱGST�����ļ�дΪg��
                </p>

                <p>��z������B�����<b>���ཻ��</b>�����õ���</b><font color=red>��=
                    90-B[1]</font>��ͬʱ�õ�<b>A��ĳྭ(�����ྭ)��</b><font color=red>�� = 90+B[0]</font>���������ྭҪת��Ϊ�Ը�����������Ȧ����ľ��ȣ���Ӧ������ྭ��ȥ�������κ���ʱ���������κ���ʱ��һ���Ƕ�������ָ������������Ȧ��Դ����ĳྭ(�Ƕ�)��
                </p>


                <p><b>��֪���������P[]���䱴�����꣺</b><br>
                    ��(1)�Ȱ�P����תΪA������ľ���<font color=red>P[0] = P[0]-��</font><br>
                    ��(2)Ȼ����ת<font color=red>P = llrConv(P,-��)</font>��<br>
                    ��(3)��Ϊ�ڳ�����ؾ��Ȳ��������ߣ����������AΪ�����㣬���������ȡ����ͨ��������ı���������ֱ�������ʽ�����Ի�Ӧ����<font
                            color=red>P = llr2xyz(P)</font>��<br>
                    ������֪���������꣬����ı������ꡣ���������輴�ɡ�
                </p>

                <p><b>��֪���屴������P[]�������꣺</b><br>
                    ��(1)P = xyz2llr(P)��<br>
                    ��(2)P = llrConv(P,��)<br>
                    ��(3)P[0] = P[0]+��<br>
                </p>

                <p><b>��ֱ�������б��������任</b></p>
                <center><img src="img/a008.gif"></center>
                <p>
                    xyz�Ǳ����ֱ�����ϵ��XYZ��A�㣨���ཻ���ϣ�����ĳ��ֱ�����ϵ��x����X���غϣ�y������(0,P,Q),z������(0,-Q,P)������<font
                        color=red>P = cos�ƣ�Q = sin��</font></p>
                <table border=1>
                    <tr>
                        <td>���(A������)ת�����</td>
                        <td>���任��Ϊ(�ⷽ�����)</td>
                        <td>�ֵ�תA��</td>
                        <td>A��ת�ֵ�</td>
                    </tr>
                    <tr>
                        <td><font color=red>x = (X,Y,Z)��(1, 0,0) = X<br>y = (X,Y,Z)��(0, P,Q) = PY + QZ<br>z =
                            (X,Y,Z)��(0,-Q,P) =-QY + PZ</font></td>
                        <td><font color=red>X = x<br>Y = Py - Qz<br>Z = Qy + Pz</font></td>
                        <td>X = X'cos�� + Y'sin��<br>Y =-X'sin�� + Y'cos��<br>Z = Z'</td>
                        <td>X'= Xcos�� - Ysin��<br>Y'= Xsin�� + Ycos��<br>Z'= Z</td>
                    </tr>
                </table>


                <p><b>��������뱴�������ת����</b></p>

                <p>
                    ������һ��ʵ����������⡣���������еĵ�P(x,y,z)���ڵ�����棬z��ྭ��γ��(��,90-��)���������κ���ʱΪg���������뾶Ϊa�뼫�뾶Ϊb����zδ֪������z��ֵ��</p>
                <img src="img/a007.gif">

                <p><b>�塢�����������������</b></p>

                <p><b>1����Բ�Ķ̰뾶��</b></p>

                <p>
                    �ӱ�������Ķ���֪�����������������ε���Ľ���Ϊ��Բ�����ĳƱ�Բ����Բ�ĳ��뾶��x�ᣬ����Ϊ�������뾶a���̰뾶��y�ᣬ��������С�أ�����������ͼ�����ǵ����һ������ͼ�����������y�����������ϣ�OD���Ǳ�Բ�İ�̾�b<sub>2</sub>���������뾶Ϊa���ϱ����뾶Ϊb��
                </p>
                <img src="img/a005.gif">

                <p>����ͼ������Բ���� (x/a)<sup>2</sup>+(y/b)<sup>2</sup>=1����x = -r*cos�ƣ�y = r*sin�ƴ���ã�<br>
                    ��<font color=red>b<sub>2</sub> = r = b/sqrt(1-c*cos<sup>2</sup>��)��ʽ��c =
                        1-b<sup>2</sup>/a<sup>2</sup>=0.00669438</font><br>
                    ����Լ���ڡ�90-̫����γ������cos�ơ�sin�ġ�sin��*sin��(���Ц�Ϊ�Ƴཻ�ǣ���Ϊ̫���ƾ�����Ϊ̫����γ)����c*cos<sup>2</sup>���Ǹ�С���������ʽ�ɸ�дΪ��<br>
                    ��<font color=red>b<sub>2</sub> = b*(1+c*sin<sup>2</sup>��/2) = b(1+0.0005296sin<sup>2</sup>��)�����е���ƫ��f=1/298.25642��b
                        = 1-f����=84381.448����</font><br>
                    ����Ӱ��պ�ɨ������Ե��Ϊ���ж����������ߣ�<font color=red>b<sub>2</sub>ͨ�����������ֵ0.9972���ɣ�</font>����ȡ��ֵ����ʵ�ϼ���½ƽ��Ҳ��2ǧ�����ҡ�
                </p>


                <p><b>2�����ཻ���뱴�ƽ��ߵļн�</b></p>

                <p>
                    Ӧ�����߶�����֤����ת��������yzƽ�档����̫�������Ǻ�Զ��������Ҫ��ö࣬����������ƶ�����Ӱ�᷽���Ӱ������΢�ģ���Ӱ��ķ�����Ҫ��̫������������̫����ͬ����ˣ����̫����γΪ��ʱ����ͼ��yӦ���ڱ�����ߡ������˵����̫����ת����γ�仯��������y���������ƶ���</p>

                <p>
                    ����̫����γС��2������̫���Ӳ�С��10���룬���Ի����ֱ�ڻƵ��棬����Ӱ�ᾭ����򣬻��������ļнǲ�����10���롣��Ȼ��ʵ�ʼ��㲻�Ƽ����ɴ�ֱ������������һ�ָ������������任��<b>���ƽ��ǣ�</b>ԼΪ90�ȡ�
                </p>

                <p>��ͼ�����������ꡢ�Ƶ������������ȵ�ʾ��ͼ��</p>

                <center><img src="img/a004.gif"></center>
                <p>ͼ�С�BOA�ɰ����·����õ����ɻƳ�任��ʽsin�� = sin�� cos�� + cos�� sin��
                    sin�ˣ�A���γ�Ǧ�=0��B��ƾ���=�ԣ������tan�� = -tan�� sin�� = -tan��
                    sin(��)����Բ������ֱ�ڻƵ�������A�Ļ�γ�½���Ϊ��BOA������<font color=red>tan��BOA = -tan�� sin��</font>
                </p>

                <p><b>3����Ӱ�뱴Բ�Ľ���(����Բ����,��Բ�õ���ƽ��뾶)</b></p>
                <center><img src="img/a010.gif"></center>
                <p><b>4����Ӱ�뱴Բ�Ľ���(����Բ����,��Բ������Բ)</b></p>

                <p>���������Բ��Բ�Ľ������ͬ����Ҫ����һ���ظ����㡣</p>

                <p><b>5����Ӱ�뱴Բ�Ľ���(��ȷ��)</b></p>

                <p><b>6����Ӱ�뱴Բ���е����</b></p>
                <center><img src="img/a011.gif"></center>

                <p><b>���ڱ�������ϵ�����������ĳ����ٶȣ�</b></p>
                <center><img src="img/a006.gif"></center>


                <p><b>��ĳһ�ر�(������)�ı�������</b></p>

                <p><b>��Ӱ����ĳ���Ӧ�ı�Ӱ�Ͱ�Ӱ�뾶</b></p>

                <p><b>Ӱ����������ڹ۲����ٶ�</b></p>

                <p><b>��ĳʱ�̵����ĵ�</b></p>

                <p><b>���ĵ�ֱ�ߵľ���</b></p>
                <center><img src="img/a015.gif"></center>
                <p><b>�������</b></p>

                <p>
                    �ڱ��������У������㷨��ֵ����¾���(���Բ�ֵ)�������һ������������ж��κ���������ͼ��Ӱ����O�˶���A��ʱT�գ����ٶ�Ϊ��</p>
                <center><img src="img/a016.gif"></center>
                <p>Rȡ����ƽ��ֵ���Ե������뾶Ϊ��λ����R=385000/6378=60.36����Ϊ���»ƾ����ƽ���ٶ�0.2128/�գ�����f =
                    1.366(T<sup>2</sup>-t<sup>2</sup>)��t�����е�O����ġ����t��B�����㣬B�˶���A��ʱT����<font color=red>f =
                        1.366t(T-t)</font></p>


                <p><b>�ճ���ûʳ���ߵĻ�ͼ����</b></p>

                <p>����ĵ��������޵ģ����Ի������ճ���ûʳ���߲���պþ����ϱ����ߵ�ͷβ�˵㣬�����㷨��������⡣</p>
                <center><img src="img/a014.gif"></center>
                <p>
                    �����ǰ�Ӱֱ������ֱ��Ӱ�ӵ�����ٶȷ���Ӱ����Խ����ϵĹ۲��ߵ��ٶȣ����������ת�йأ���Ӱ���ٶȱȱ����ƽ���ϵ�����ת�ٶȴ���࣬���Ժ���ʵ�ʻ�ֱ��Ӱ�ٷ�����б���ࡣ����Ӱ�ӵ��ƶ��������뱴Բ������Ľ��㹹���ճ�ûʳ���ߡ�</p>

                <p>
                    ��ͼ�е�����ͼ�������粻����ʱ��q1����д�����q3�����ڣ��������ǰ����һ��q3������q3��ͷд��q1�����һ�㣻���������ʱ��q1��β��q3��ͷ�ֱ��ϱ������ͷ����β�������Ч������������һȷ��q1��q3��Ϊ���ߣ����ʹ��ͷ���ӵ�����Ķ˵㡣q2��q4����ķ�����ͬ��</p>

                <p><b>��ͼ������</b>ֻҪ������ڣ�����q1��q3�Ƿ���ڣ�q1��βǿ�ƽӱ���ͷ��q3��ͷǿ�ƽӱ���β�������ڣ�q3�ز����ڣ���ǿ�Ʋ�һ����q1�����������һ��q3���Ͻ�Ĵ��?�����ơ����q1��q2��ͷ���ӣ�q3��q4��β���ӡ�<br>
                    &nbsp;&nbsp;1��q1��q2��q3��q4�ֱ�洢����ʳ���ߣ�q1��q3Ϊʳ�����ߣ�q2��q4Ϊʳ�����ߡ�����дָ��Aָ��q1��q3������дָ��Bָ��q2��q4����ʼ״̬A=q1��B=q2<br>
                    &nbsp;&nbsp;2��д�����У���ʱ����Խ�����ʳʱ�̣� ���q1��������A�Ƶ�q3�����q2��������B�Ƶ�q4��<br>
                    &nbsp;&nbsp;3��д�����У����б����ж��ߣ��޽⣩�����A�Ƶ�q3��������߶��ߣ���B�Ƶ�q4��<br>
                    &nbsp;&nbsp;4�����߶�q1��q2��q3��q4ȫ��д��֮��(1)����ޱ�����q3��Ϊ�գ���q3����һ��q1��ĩ��꣬ʹq1��q3���������������ʳ����q1��q3�����Ϊ�գ�q1ĩ����һ������ͷ��꣬q3������һ������β��ꡣ����q1��q3��������ı�ʳ���ߡ�(2)q2��q4��ͬ���?ʹ��q2��q4����������ʳ����<br>
                    &nbsp;&nbsp;5���ϡ���ʳ���ߵ���ͷ����β�ֱ𣬼�q1��q2��ͷ���ӣ�q3��q4��β���ӡ�
                </p>

            </td>
        </tr>
    </table>
</center>
</body>
</html>
